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Abstract 

We consider time-harmonic wave scattering from an inhomogeneous 
isotropic medium supported in a bounded domain n cR^ {N > 2). In 
a subregion L) d fi, the medium is supposed to be lossy and have a large 
mass density. We study the asymptotic development of the wave field as 
the mass density p — )■ +00 and show that the wave field inside D will 
decay exponentially while the wave filed outside the medium will converge 
to the one corresponding to a sound-hard obstacle D <£ Q buried in the 
medium supported in 0,\D. Moreover, the normal velocity of the wave 
field on dD from outside D is shown to be vanishing as p — )• -|-oo. We 
derive very accurate estimates for the wave field inside and outside D and 
on dD in terms of p, and show that the asymptotic estimates are sharp. 
The implication of the obtained results is given for an inverse scattering 
problem of reconstructing a complex scatterer. 



*Departnicnt of Mathematics and Statistics, University of North Carolina, Charlotte, NC 
28263, USA. (hliu28@uncc.edu) 

^Institute of Mathematics, Academy of Mathematics and Systems Science, Chinese 
Academy of Sciences, Beijing 100190, P. R. The work of this author was partially supported 
by grant under NSF No. 10990012. (zaijiu@amss.ac.cn) 

■^Institute of Mathematics, Academy of Mathematics and Systems Science, Chinese 
Academy of Sciences, Beijing 100190, P. R. China. The work of this author was partially 
supported by grant under NSF No. 10990012. (hpsun@ainss.ac.cn) 

^Department of Mathematics, Chinese University of Hong Kong, Shatin, N. T., Hong 
Kong. The work of this author was substantially supported by Hong Kong RGC grants 
(Projects 405110 and 404611). (zou@math.cuhk.edu.hk) 



1 



Singular perturbation of the reduced wave equation 



2 



1 Introduction 

We shall be concerned in this paper with the following scalar wave equation 
(see, e.g., 



c^(x) dt"^ dt \p{x 

for all X eR^ {N >2) and t eR+. In equation (JTI]), U{x, t) is the wave field, 
c{x), (j{x) and p{x) are positive scalar functions and represent the wave velocity, 
the damping coefficient and the mass density of the medium respectively. It is 
supposed that the medium is compactly supported in a bounded domain Q in 
M^. We consider the medium outside Q to be homogeneous and no damping 
present, so we may assume after normalization that c = cq, p = 1 and a = 
in Q'^ := M^\n. Let D (£ Q he a subregion of Q and the material parameters 
inside D be given by 

c(x) = Co, cr{x) = (To, p{x) = for x E D , (1.2) 

where Co,o"o and e are positive constants. This work shall be devoted to 
the study of the asymptotic development of the wave field U{x, t) as the mass 
density p inside D tends to infinity, i.e., the parameter e ^ 0~^. We shall consider 
the time-harmonic wave propagation, namely to seek a solution of (11.11) in the 
following form 

U{x, t) = 3?{M(x)e~^^*}, F(x, t) = 3?{/(x)e"*'^*}, 

where u E M+ is the frequency. By our earlier assumption on the homogeneous 
space outside the medium Q, we see the wave number k = u/cq. We suppose 
that f{x) is compactly supported outside the inhomogeneous medium, namely 
supp{f) C Br^\Q for some i?o > 0, where and in the sequel Br denotes a ball 
of radius r centered at the origin in M^. Factorizing out the time-dependent 
part, the wave equation (11. ip reduces to the following time-harmonic equation: 

V . ( -Wu] +k'(^ + t^] u = fix) in M^. (1.3) 



We shall seek the total wave field of (II. 3p admitting the following asymptotic 
development as |x| — )■ oo: 

u{x) = e^'-' + ^^^^ Ia (x, d,k) + I , (1.4) 
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where e^^^''^ is the incident field, and A{x, d, k) with x = x/\x\ is known as the 
scattering amphtude (cf. [2] [5]), with d G For notational convenience, 

we set 

and u'^(x) = — is the scattered field outside the medium region Q. 
Throughout the rest of the paper, we assume that fl and D are both bounded 
domains such that R^\U and n\D are connected. Let q e L°°{Q\D) and 
7(x) G satisfying the following physically meaningful conditions: 

70 < 7(a;) < To, 3f?g(x) > To, '^q{x) > ioi x e n\D , 

where 70, Tq, Fq are positive constants. With all these preparations, we can 
formulate our interested problem of finding the total wave field u{x) of form 
fOj) to the system ([L3D as follows: Find Ue G H^M^) such that 



V ■ (sVue) + k\r]o + iTo)Ue = 




in 




V ■ (7(x)VMe) + k'^q{x)ue = 




in 


n\D, 


An^ + k\'^ = f 




in 




Ue = u"- + ul 

, du" du't 
, ,• 9m, dut 

hm |a:|(^"i)/2|f''', z^l 


dv 
= 0, 


in 

on 

on 


91], 



(1.5) 



where u denotes the exterior unit normal to dD or dQ. We use the notations 
uj to represent the limits oiu^ on dD or dfl, taking respectively from inside 
and outside D or Q. The last limit in (11. 5p is known as the Sommerfeld radiation 
condition. The well-posedness of the scattering problem (II. 5p is given in the 
Appendix and the scattering amplitude in (II. 4p can be read off from the large 
asymptotics of u^. It is readily seen that Us depends on e nonlinearly and so 
does u^. In order to present the main results of this paper, we introduce the 
following scattering problem: 
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Find u e Hl^(R^\D) such that 



V ■ (7(x)Vm) + k'^q{x)u = 
An" + 



f 



u = u 

u~ 



u 



u + u\ 7 



hm 



X 



(N-l)/2 



du 
du 

d\x\ 



du . du 



in n\D, 
in M^\n, 
in M^\n, 

on dD, 
on dQ, 



;i.6) 



— iku' 



0. 



One can see from ( \1.6\} that the normal velocity of the wave field vanishes 
on the boundary dD, so the wave can not penetrate inside D. In the acoustic 
scattering, D is known as a sound-hard obstacle, so the system f ll.6p is an 
obstacle scattering problem with an obstacle buried inside some inhomogeneous 
medium. We shall show that the solution of the medium scattering problem 
(11. 5p will converge to the solution u of the obstacle scattering problem (11. 6p as 
£ ^ 0+, or the density p of the medium D tends to infinity. This is reflected by 
the results in the following three theorems, where C and C are generic constants, 
which depend only on q, k, rjQ, tq, 7, Sq, D, fl, Br, but completely independent of 
e. 

Theorem 1.1. LetUe E if/„^(M^) andu E Hl^{R^\D) be the solutions to / fO) 
and U.6\) . respectively. Then for any R > Rq, there exist Eq > and C > 
such that the following estimate holds for e < Eq: 



F \\m{BR\n) + 



LHBR^\n))- (1-7) 

As a consequence, the scattering amplitude As of converges to the amplitude 
A of in the following sense that 



for some constant C > and all e < Eq. 



\m{Bii\n) 



L^Bu^\n)j 



The next theorem characterizes the normal velocity of the wave field on 
the boundary of the medium D. 

Theorem 1.2. For the solution E Hl^^iM'^) to the system ( fi.5]) . there exists 
£0 > such that the following estimate holds for e < Eq: 

""'""miBnSn) ^ \\J \\L^BRg\n)j 



7- 



du 



;i.9) 
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Moreover, the next lemma indicates that the solution inside the medium 
D decays exponentially. 

Theorem 1.3. Let Dq be a suhdomain such that Dq D with dist{dDQ,dD) 
> So > 0, and yjrjo + ztq = a + bi with a > 0,6 > 0. Then for the solution 
G Hl^J^^) to the system M.5\] . there exists £0 > such that for e < Eq, 

khS 

\\ue\\c(Do) < C'exp(-— -=) i\\u'\\HHBR\n) + \\f\\L^BR \n) 

). (1.10) 



2 Discussions 



We are interested in the scattering from a compactly supported inhomogeneous 
isotropic medium, with a subregion occupied by some medium possessing a large 
density. Based on our discussions in the previous section, we let 

{n\D;^,q}®{D;e,rio + iTo} (2.1) 

denote the inhomogeneity supported in Q in (11. 5p . and 

{n\D;^,q}(BD (2.2) 

denote the scatterer in (II. 6p . where D is known as an impenetrable sound-hard 
obstacle in the acoustic scattering (cf. [2]). As it can be seen from (11.61) . the 
wave field for a sound-hard obstacle can not penetrate inside and the normal 
wave velocity vanishes on the exterior boundary of the obstacle. We call the 
scatterer in (12. 2p . composed of an obstacle and a surrounding inhomogeneous 
medium as a complex scatterer. In this work, we actually show that 

{n\D;j,q}®{D;e,r]o + iTo} {n\D;-f,q}®D as e ^ 0+, (2.3) 

in the sense of Theorems I1.1H1.3[ That is, a sound-hard obstacle can be treated 
as a medium with extreme material property, namely with a very large mass den- 
sity. Despite the nonlinear nature of the convergence (12. Sp . we can still derive 
very accurate estimates in a general setting. In addition to provide a mathe- 
matical characterization of a physically sound-hard obstacle and its asymptotic 
connection to media with extreme material properties, we would like to note 
that the results established in this work could have some interesting implication 
in the inverse scattering problem of reconstructing a complex scatterer. In fact, 
it can be seen that a complex scatterer could be reconstructed as a medium, and 
one could locate the embedded obstacle in the reconstruction as the subregion 
with a large density parameter. 
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Finally, we make another practically meaningful remark on our study. In 
(12. ip . the outer inhomogeneous medium {Q\D]'~f,q} could be anisotropic, for 
which one could also show the convergence ( 12.3P by modifying our arguments 
in the subsequent sections. However, as mentioned earlier, one of our main 
motivations is from the inverse scattering problem. If the surrounding medium 
is anisotropic, one could not uniquely recover a complex scatterer; actually 
one may have the invisibility or virtual reshaping phenomena (see, e.g. [8] 
[1]). This is why we focus on the isotropic setting in this work. The extreme 
medium inside D is assumed to be lossy, which is a realistic assumption from 
the practical viewpoint. 

The rest of the paper is organized as follows. In Section [3|, we prove 
the main results of this work, and demonstrate the sharpness of our major 
theoretical estimates by considering a special case based on series expansions in 
Section HI 



3 Proofs of the main theorems 



This section is devoted to the proofs of Theorems I1.1H1.3I in Section [H For the 
purpose we need the following lemma. 

Lemma 3.1. Consider the following transmission problem 



V ■ {'y{x)Vv) + k'^q{x)v = 
Au' + k\' = f 

du 



in Q\D, 



in 



du du 
hm 



on dD, 
on dQ, 
on dfl, 



{3.1] 



d\x\ 



iku'' 



0. 



There exists a unique solution {v,u'^) G H^{Q\D) x Hl-^^(R'^\Q) to liS. and 
the solution satisfies 



\\''^\\m{n\D) + W'^^^Wh^BrXq) 
<C'(lblli/-i/2(0D) + llfi'i||/fi/2(an) + llfi'2||/f-i/2(an) + \\j \\L^{BR^^\n)) 



(3.2) 



where the positive constant C depends only on 7, q, k, Q, D and Br, hut inde- 
pendent ofp,gi,g2,f. 
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We could not find some references on tlie well-posedness of the transmission 
problem (13.11) , so provide a proof by using a variational technique presented in 
[3] and [H]. We first demonstrate the following auxiliary lemma. 

Lemma 3.2. The system fl3.ip is uniquely solvable and it is equivalent to the 
following truncated system: find {vi,Ui) G H^{Q\D) x H^(Bji\Q) such that 



V ■ (7(x)Vfi) + k^q{x)vi 
Aui + k^Ui = f 
P 

vi-ui = gi 

92 







I du 



' dv du 



in Q\D, 

in BR\n, 

on dD, 

on dQ, 

on dQ, 

on OBr, 



(3.3) 



where A : H^^'^(dBji) — H ^/"^{dBji) is the Dirichlet-to-Neumann map defined 
by A^jj = f^las^ (cfi /f /f with W e HI^^{R^\Br) bemg the umque 

solution to the system 



AW + eW = 

W = i/je H^/^idBn) 



in R^\Br, 
on dBn, 



lim 



\x\- 



\X 



(7V-l)/2 



dW 
d\x\ 



ikW 



(3.4) 



0. 



Proof. We first show the uniqueness of the solution {v,u'^) to system fl3.ll) . For 
the purpose we set p, gi, g2, f to be all zeros. Multiplying the first and second 
equations of (13. ip . respectively, by v and u^, and integrating by parts in Q\D 
and Br\Q, together with the use of the boundary conditions on dD and dQ, 
we have 



'j\Vv\'^dx + 



n\D 



k'^q\v\^dx 



n\D 



Wu'l'^dx 



BR\n 



+ 



BR\n 



du 



(3.5) 



u'^ds = 0. 



Taking the imaginary part of both sides of (13. 5p . we derive 



du 



■u^ds 



k q\v\ dx < 0. 



n\D 



Then by Rellich's lemma (cf. [2]) we know is zero outside Br, which with 
the unique continuation implies that = in Q\D and v = in D. 
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Next we show the equivalence between systems (13.11) and (13.31) . By the 
definition of A, we see that if {v,u'^) solves the system (13. ip . then (fi = v,ui = 
u^Ibj^Xu) is the solution to the system ( 13. 3p . On the other hand, by applying the 
Green's representation (cf. |2](2.4)) to the solution {vi,ui) of (13. 3 p we obtain 
that 

+ / (Au,iy)^x,y)-uM^^^) ds{y) - [ f{yMx,y)dy, 



for X E Bfj\Q, where 



^ X (7V-2)/2 



= 4 ^^^TF^J HlZ,y,ik\^-y\) (3.7) 

is the outgoing Green's function. By definition of A and the radiation of $(a;, y) 
(cf. pp. 98 in [3], and [H]) 



AmivMx, y) - ni(y)^|^ ) ds{y) = 0. 



Hence, 



u,{x) = -f f^^<f(x,y)-ni(y)^|^')ds(y)- / f{y)^{x,y)dy. 

(3.8) 

It is clear that ui can be readily extended to an Hl^J^^\VL) function, which 
we still denote by ui. We can see that ui satisfies the Sommerfeld radiation 
condition, which together with the uniqueness of solution to (13. ip implies that 

Ml = U'^. □ 



With the uniqueness and equivalence in Lemma 13.21 we can apply the vari- 
ational technique to study the reduced problem (13. 3p to prove Lemma 13.11 

Proof of Lemma \3.1\ Without of loss generality, we assume fc^ is not a Dirichlet 
eigenvalue of —A in B{i\Q, and introduce the following auxiliary system 

-Af; -k^ = in 
V = gi on Sfi, (3.9) 

f; = on dBfj. 
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It is easy to see ||^'||Hi(Bjj\fI) — ^llfl'illHi/2(af7). We now set 

jvi{x), xen\D, 
w{x) := s _ — 

I + v{x), X e Bfi\Q. 

We can check that w G H^{Bfi) satisfies the following equation: 



V ■ (7(x)Vm;) + k'^q{x)w = 
Aw + k'^w = f 

7f =P 



w = 

' du 

dw 

< dv ~ 



du 



dv 

92- 



av 



in Vl\D, 
in BR\n, 
on dD^ 
on dVt, 
on dVt, 
on dB^. 



Next, we define Aq: H^/^dBR) ^ H-^/'^{dBR) by 



AoV'i 



du 



dBf 



where Wi G HI^^{W^\Br) is the unique solution of the system: 



-AWi = 



in R^\Br, 
on OBr, 



(3.10) 



(3.1i: 



and satisfies the decay property at infinity, namely Wi = 0{\x\~^) for 
and Wi = C(log for iV = 2, as |a;| — t- +oo. 
It is known that (cf. [3] and [TT] ) 



(3.12) 

iV = 3, 



-/ Mo^iC?s>0, V^i G i/^/'(95^), (3.13) 

and A — Aq is compact from H^^^(dBR) to H~^^^{dBR). Then for any (/) G 
H^{Br), using the test function (f we can easily derive the variational formula- 
tion of system (13. lip : find w G H^{Br) such that 



ai{w, ip) + a2{w, if) = J^{ip) 



(3.14) 
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where the bihnear forms ai and 02 and the hnear functional J-" are given by 

ai{w,(p):= / 7VU' ■ V(^(iy + / k'^wcpdy + / Vw-Vipdy 
Jn\D Jn\D JbrXq 

k^WLpdy — I AoWifds, (3.15) 

Bii\n JdBn 

a2{w, if) := — / k'^{q + l)w(pdy — 2 / k'^w^pdy — I (A — AQ)w(pds, 

Jn\D JbrSQ JdBR 

(3.16) 

^(V^) := - / P'^ds + (5(2 - ■^)'^ds + / -^vds - / f<^dy. 
JdD Jan du Jqb^ du Jb^ 

(3.17) 

Using fl3.13p we can readily verify that for any (pjip E H^{Br), 

\ai{(t),^)\<Ci\\(p\\H\BR)\W\\m(BR) and ai((p, ^) > C2||(^||^i(^^) (3.18) 

for some constants Ci and C2. Then by Lax-Milgram lemma there exists a 
bounded operator £ : H^{Bji) — )• H^{Bji) such that 

a,{w,ip) = {Cw,^), W^,weH\Bj,), (3.19) 

where and in the following, (-, ■) denotes the inner product in H^{Br). Moreover, 
the inverse exists and is bounded. By Riesz representation theorem, we also 
know that there exist bounded operators /Ci, /C2 : H^{Bf>) — )■ H^{Bfj) such that 

03(1^,^9):= / k'^{q + l)wLpdy + 2 / k'^wipdy = {JCiW, cp) (3.20) 
Jn\D JbrXu 

and 

ai{w,(p):= / {A — Ao)wipds = {)C2W, (f). (3.21) 

JdBR 

We now claim that both /Ci and /C2 are compact. In fact, let {w„}„6N be a 
bounded sequence in H^{Bji) and < M, and we can assume that 

Wn Wo in H^{Bjij. Since H^{Bji) L'^{Br) is compact, we know Wn Wq 
in L'^{Bn). By we can write 

a-siwn - Wo, if) = {lCi{wn - Wo), (p). (3.22) 

Taking ip = /Ci(w„ — Wq) and using (13.201) . we can verify that 

\\ICi{Wn-Wo)\\m{BR) < 4:Mk'^max{\\\q+l\\j^^^^^\-^^,2}\\}Ci\\\\Wn-Wo\\L2{BR) ^ 0, 



Singular perturbation of the reduced wave equation 



11 



which imphes the compactness of /Ci. In a similar manner, we can prove the 
compactness of /C2. Indeed, let Wn wo in H^{Br), and by trace theorem, 
WnldB^ ^ wolaBn in H^'^dBn). Since A - Aq : H^'^dBn) ^ R-'/^dBn) is 
compact, we see (A — Ao)w„ — > (A — Ao)wo in H~^^'^{dBji). By ( 13.2 ip we can 
write 

Taking ip = JC2{wn — wq) and using (13.211) . one has 

\\lC2{Wn - Wo)\\hHBr) < II(A - Ao)K^ - Wo)\\H-i/2(^dBn)\\^2{Wn - W^o) ||iTi/2(aiJfl) 

< CsMlKA - Ao){Wn - Wo)\\H-^/^idBn)\\J^2\\ ^ 0, 

which implies the compactness of IC2- 

Since £ is bounded and invertible, and /Ci + IC2 is compact, we know C — 
(/Ci + JC2) is a Fredholm operator of index zero. By the uniqueness of (13. ip . 
(£ — (/Ci + /C2))~^ is bounded. On the other hand, it is straightforward to show 

\F{ip)\ < C{\\p\\fj-i/2(^Qj;)) + \\gi\\Hi/2^QQ) + \\g2\\H'y^dn)Hf\\L^BR^+^\BR^)^^^^ 
which readily implies (13. 2p . □ 

The next lemma presents some important a priori estimates of the solution 
Us to (II. 5p in terms of e. 

Lemma 3.3. Let G Hl^J^^) he the unique solution to M.5\) . There exists 
£0 > such that the following estimates hold for all e < Eq, 

\\'^e\\m(BR\D) ^ C'l(ll/llL2(Bfly\?I) + lk1l/fi(Bfl\T!)) ; (3.23) 
V^heWmiD) < C2{\\f\\L2^Bn,\n) + M\HHBn\n)) (3-24) 

where the constants Ci and C2 are independent of e. 

Proof. Multiplying to the both sides of the first and second equations of (11.51) 
and integrating over f2, we have 

- / e\Vue\^dy+ I k^{rio + iTo)\ue\'^dy - / _'y\Vus\'^dy 

k^q\us\'^dy + / 7— -^m^c/s = 0. 
Q\D Jan c*// 

Then multiplying to the both sides of the third equation of (11.51) and inte- 
grating over Br\Q, we obtain 

^Kds+ [ ^Kds- [ iVu^l'dy 
an diy Jqb^ du Jb^\^ 



k'^\ul\'^dy = / /u^c 

Ba\n JBa\n 



(3.26) 
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By adding up fl3.25p and fl3.26p . using the corresponding transmission conditions 
and then taking the imaginary and real parts of the resulting equation, we derive 



k'^To\ue\'^dy + / k'^'^q\ue\'^dy + Q / —^u^ds + Q / ^:—utds 
D Jn\D Jan du Jg^ du 

$5 [ ^u'ds + 53 /" ^utds = 53 /" futdy 
Jan du Jqbj, du J^^^xn ^ 



(3.27) 



and 



- / e\Vu^\ dy+ k r/okel dy - _'j\Vue\ dy 
Jd Jd Jn\D 

+ / k^^q\uA''dy + ^ / ^u'ds + ^ / ^uids 
Jn\D Jan Jq^ du 

+ ^ f ^u'ds + ^ [ ^uids- [ iVu'fdy 
Jan ov Jqb^ du Jb^^xu 

+ / _k^\u'fdy = ^ [ _fuldy. 
From (I3.27p . one has by direct verification that 

+ \\f\\L'^{Bit\n)\\K\\m{BR\n) 

< 8C {\\Ue\\%^BR\D) + \\^'\\miBR\n) \\J \\L^Bii\n)l ' 



(3.28) 



(3.29) 



where C depends only on ?7o, tq, k, q, fl, Br. We can readily check by ( I3.28P that 

j^e\Vu,\'^dy < C2 (^|he|li2(D) + he\?Hi(BR\D) + II^II^Hi^fl 



+ \\i\\L^{BR\n)\We\\m{BR\n) I' 



\n) 

(3.30) 



where C2 depends only on k,rjQ,q,^,fl, Br. Combining f l3.29p and (13. 30 p . we 
see that there exists a constant C3 dependent only on k, q, 779, Tq, 7, fi, Br, such 
that for e < 1, 

— / \ 1/2 
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Next, we prove (13.231) by contradiction. Suppose fl3.23p is not true, then 
without loss of generahty, we can assume that for each n G N, there exist 



and Un such that ||/'"||x,2(5 



n\\H'^{Bii\n) = 1 the corresponding 



solution tends to infinity, i.e., ~^ +^ as e — )■ 0+. Let 



u: 



u 



\u 



fe 



e \\m{Bii\D) 

r 

\\m{BR\D) 



\\^e\\m{Ba\D) 
W'e\\m{BR\D) 



(3.32) 



Clearly, Ve,n e H^, 

and the source We have 



is the unique solution of (ll.Sp with the incident wave 



(3.33) 



By a completely similar argument as we did in deriving (I3.3ip . we can show 
that for sufficiently large n, 



s,n\\H^D) 



<C, 111;, 



^e,n\\Hi^BR\D) 



i ||2 



e,nnH^Ba\n) 



+ Wfe IIl2(Bh\!^) 



1/2 



<C.V2. 



(3.34) 



By taking the trace and using the transmission condition on dD and (I3.34p . we 
know the existence of a constant C4 depending only on D such that 



dv: 



7- 



du 



dv: 



dv 



(3.35) 



Noting that (fe,n 
91 = K.nlan, 92 = 



n\D' '^e,n\ 
di'i.n I 



dv 



is the unique solution of (13. ip with p = 7 
then by Lemma 13.11 we have 



dvj.n 

du 



\V. 



s,n\\H^{Bji\D) 



< c 



7- 



dv 



+ 11/; 



e \\L^{Bn\n) 



By (13331) . (13351) and (l336l) . we further derive 



(3.36) 



which contradicts with the equality \\ve,n\\H^(^BR\D) — ^ ^'^'^ ^^^^ proves (I3.23p . 
Now by combining (l3:23|) with disT]), we obtain (E^ll)- □ 
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We are now in a position to present the proofs of Theorems ll.H - ll.3[ 



Proof of Theorem This is a direct consequence of (I3.24p in Lemma 13.31 
Indeed, by taking the trace on dD, we see 



duz 



du 



< C\\uA\m 



s\\m{D), 



where C depends only on D. Then by the transmission condition on dD, we 
readily derive ( II .Qp : 



7- 



du 



f/-i/2{aD) 



du 



+ W 



\m{Bii\n) 



□ 



Proof of Theorem \l.l\ Let V = — u, V"^ = ul — . One can verify directly 
that V satisfies equation (13. ip with f = 0, p = 7^ = 7^^|az) and gi = g2 = 0. 
Then by Lemma 13.11 and Theorem 11.21 we have 



<Cey' {\\f\\LHB,,\n) + hlHHBn\^)- 
Finally we know from [2] (pp.21) that 

r (^,e~'^^-y dv 

'dBR I 



7- 



dv 



H-i/2(aD) 



(3.37) 



{As - A) (x) = C 



ikx-y 



dv 



dv 



ds{y), X E 



(3.38) 



where ( = I/Att for = 3 and C, 
one can derive (II. 8p by some straightforward estimates. 



; for iV = 2. Using (l337j) and (l338|l . 



□ 



Proof of Theorem \1.3[ We shall make use of the following integral representa- 
tion of the wave field inside D (cf. [2]): 



uAx 



dD 



ds{y), X E D, (3.39) 



where G{x, y) is the fundamental solution corresponding to the first equation of 
(II. 5p and is given by 



G{x,y) 



^ik\x-y\ 

47r|x — y\ 



for iV = 3; G{x,y) = '^Hi^\k\x - y\) for iV = 2 , 



(3.40) 
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with k = k{a + ib)e~^^'^. 

Next, we shall only prove the theorem for the 3D case and the 2D case could 
be proved in a similar manner. For x G -Do and y G dD, since |x — 1/| > 60, it 
can be verified by straightforward calculations that 



Jk\x-y\ 



4:7l\x — y\ 
^ik\x-y\ 



47r|x — y\ 



< 



< 



47r(5n 



kVa^ + 62 
^172 



1 

To 



On the other hand, by fl3.24p in Lemma 13.31 we see that 



du 



\m/^(dD) 



H-i/2(az)) 



{3A1] 



(3.42) 



Now using fl3.4ip and f l3.42p in f l3.39p . one can obtain fll.lOp by straightforward 
calculations. □ 



A special case and sharpness of convergence 
estimates 



In this section, we shall consider a special case of the model system (11.51) : D 
is the ball B^i^ of radius Ri, and only the subregion D is occupied by the 
inhomogeneous medium in the whole space , and the rest is the homogeneous 
background, so we have 7 = 1 and g = 1 in (II. 5p . Moreover, we consider the 
scattering only from plane wave incidence, namely, / = 0. We shall derive the 
corresponding estimates of the wave field, which shall demonstrate the sharpness 
of our convergence estimates in Section [3l We will consider only the 3D case 
while the 2D case could be treated in a similar manner. 

In our current special setting, we can rewrite the equation (II. 5p as follows: 
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Find Us{x) G -^/^^(M^) which solves the system 



V ■ (eVue) + k'^ivo + iTo)ue = in D, 



Auj + k^Ue = 



^ikx-d _|_ u^^l^x) 



lim \x\ 

\x\—¥oo 



in M?\D, 
in M?\D, 

on 



(4.1^ 



IS - 1^ = 0, 

o\x\ 



and the equation fll.6p with D a.s a. sound-hard obstacle reduces to 
'Au + = 



lim I X I 



in M3\L), 
in M3\D, 

on dD, 



(4.2) 



- iku" 



In the sequel, we let go = (^o + i^o)/^ and y/qo = e~^^'^{a + bi) with a > 0, 
b > 0. We shall make use of the spherical wave series expansions of the wave 
fields in (14.11) and (14. 2p . and we refer to [2] for a detailed discussion about 
spherical wave functions. Let ^^(x) and m| be given by the following series: 



oo n 



U^iX] 



n=0 m=—n 
oo n 



(4.3) 



n=0 m=—n 

where x = x/\x\, and u'*(x) be given by 

oo n 
n=0 m=—n 

We shall make use of the following series representation of the plane wave 



(4.4) 



oo n 



ikx- 



E ^"4vrF-(rf)j„(A:|x|)F; 



(4.5) 



n=0 m=—n 
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By (14.31) and (14.51) . and using the boundary condition on dD, we know 



m 



Next, by the transmission boundary conditions in (14.11) and comparing the 
coefficients of YI^(x) we derive 



ek^oh'^mV^oRi) = ka^h^^y{kR,) + t^k47rY^^id)f^{kR,). 
Solving the equation (14. 6p . we obtain 



t^A7rY;^{d)j'^{kR,)Uk^R,) - e^t^47TY;^{d)f^{k^R,)jn{kR, 



e^j'n{k^Ri)h'k\kRi) - h'^^y {kR^)3n{k^oR 



-i''A'KY^{d)jn{kRi)h^n'^\kR^) + i'^A7TY;^{d)h'n\kRi)f^{kR 
e^o3n{k,/qoRiWn\kRi) - h^n'^' {kRi)jn{k^Ri) 



(4.7) 



We first consider two wave fields outside D and show the following lemma, 
which indicates the sharpness of the estimates in Theorem 11.11 

Lemma 4.1. For the far field patterns As and A corresponding to the solutions 
Us and u of systems OTTl) and we have 



\Aeix)-Aix)\=CAe^^^ + 0{e), Vx G (4.^ 
where C_a depends only on rio, tq, k, Ri, d. 
Proof. In fact, by (14. 3 p and (14. 4p . we have 



n=0 m=—n 
oo 



(4.9) 



But it follows from ( K7\\ that 



_ z^A7rY^{d)UkR,) - T{qo,n)z-A7rY^^{d)jn{kR,) 
T{qo,n)hL'\kRi)-hl^^ (kR^) 
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with 



T{qo,n) := 



jn{k^Ri)' 

Next, we derive the asymptotic development of T(go, n) as e — ?■ 0+. Noting that 
fni^) = -jn{z) -Jn+i{z) (cf. [2]), we See 



j'nikVqoRi) 

then 

T(go,n) =ey/% 



n 



n 



ky/%Ri 

jn+i{ky/q^Ri 



3n{ky/q^Rl) - jn+l{ky/%Rl] 



ky/q^Ri jn{ky/%Ri] 



ne ^jn+iik^/q^Ri 



jn{ky/%Rl) 

(4.11) 



In virtue of the asymptotic behavior of jn{z) (cf. 9.2.1 and 10.1.1 jT]) as 

\z\ — > oo and lare; z\ < vr, one has 



j^{z) = -{cos(2 - mT/2 - 7r/2) + e\'^'\0{\z\-')} 

and as e — J- +0 (cf. [^), one also has 

jn+iiky/qi^Ri) 



J7T/2 



jn{ky/%Rl) 

Combining ( I4.1ip - p.l3p . one has by direct calculations 

Jn+i{ky/q^Ri 



(4.12) 



(4.13) 





ne 




< 




+ 




kRi 



\nqo,n)\ < 



Now, by (14. 9p . we have 



jn{ky/%Ri] 



0{ne + y/I). (4.14) 



oo n 



-1 



?i=0 m=—n 



In the sequel, we let 



h^n^ {kR,) 

j'nikRl) 



h^n^\kRi) 



By using the Wronskian jn{t)y'^{t) — f^{t)yn(t) = 1/t^, we then have 



zT(go,ri)47ri;-(c/) 



k^Rl[T{qo,n)hl^\kR,) - h^^^ {kR,)]h^^^ {kR,) 
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Next by the asymptotic behavior of hn\kRi) (cf. [2]), 



1 ■3---(2n- 1), ^^,1, 

n 



and also using the relation hn^ [z] = —h^n+ii^) + —hn\z), we have 



47rF™(c?) 



(4.16) 



By (14.161) and (14.181) . one readily sees that for sufficiently large n and small 



K^RlKn^ {kR^Y 
so constant in (14.81) can be chosen as 



(4.17) 



oo n 



n=l m=—n 



Noting that for any n, m G N (cf. [2]) 



|y-(d)ir(i:)l < 



2ra + 1 
47r ' 



(4.18) 



hence is bounded. Finally, using (I4.17P and the asymptotic development 

of hn^ (kRi) for large n (cf. [2]), one can show (14. 8 P from (14.150 by direct 
calculations. □ 

Next, we consider the normal velocity of the wave field on and show 
that there exists a constant C^, which depends only on k, Ri, d, rjo, tq such that 



duZ 



du 



(4.19) 



Clearly the estimate (I4.19P shows the sharpness of the estimate in Theorem 11.21 
In fact, by the transmission condition on we have 



oo n 



' -e^\dBn,=ek^oYl E KJnikV^oRi)Y: 



du dv 



[X . 



n=0 m=—n 
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Using the Wronskian relation, in{t)y'^{t) — j'^{t)ynit) = we get 



By direct calculations we obtain 



(4.20) 



vn=0 m=—n 



, , .XN -1/2 



1/2 



(4.21) 



Then by dUlD, fO0|) and the asymptotic behaviors of h'n\kRi) and /il^^ (A;i?i) 
for large n (cf. [2]), one can show that the series involved in fl4.2ip converges to 



n(n + 1' 



-1/2 



n=0 m=— n ^ 

as £ — 0^. Hence, for e sufficiently small we have 



167r^|F„'"(rf)|^ 
k^Rl\h^rl^\kR^)\^ 



duf 






dv 




du 



= C,^e + 0{e), (4.22) 

with a = 2ky%{a^ + 62)1/2. 

Finally, we consider the wave field Ue inside B^^ ^ ^Ri with So = Ri — 
i?2 > 0. By (14. 3p . it suffices for us to consider the asymptotic development of 
bnjn{ky/q()\x\) for |x| < i?2- We first note that 

Cjn(fcygokl) = b'^jnik^Ri 



^ jn{ky/%\x\) 

jn{k^Ri) 



■n+1 



47rr™(c/) 



jn{ky/%\x\) 



(4.23) 



k^Rl{T{qo,n)h^n\kRi) - h^l^' (kRi)} 3n{k^Ri 
By ( I4.12P one sees that 

|jn(fcA/goi?i)| ~ - 



Ri 



as e 0+ . 



(4.24) 



In the sequel, we consider two separate cases for Us{x) with x G Bn^- First for 
the case that |A;y^||x| = A;£:~"'"/2|a + i6| > 1, then l/\x\ < ke~^^'^\a + ib\, and 
we can show 

jn{k^\x\) 



jn{ky/%Ri 



Ri 



\x\ 



(4.25) 
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as e — )■ 0"*". Hence by combining fl4.18p . f l4.23p with f l4.25p we derive that 

oo n 

<E E Kjnik^0\x\)Y:^{x)\ 



n=0 m=—n 

oo n 
n=0 m=—n 



k^Rihil^'ikRi] 



<Mik\a + ih\e-^^^^/^^' , Vx G Br^ 
for sufficiently small e such that e~^^'^\a + i6| exp(— A;6(5o/(2a/£)) < 1, where 



n=0 m=— n 



2(2n + 1) 



k^Rlh^n'^\kRi] 



< +CX). 



For the other case, if I /c^go I = ^/^|a+i6||x| < 1, then using the asymptotic 
behavior of jn{z) for large n we know there exists a constant M2 such that 



|jn(A;v^|a:|)| < M2, VriGN. 



(4.26) 



In a similar manner as we did above one can obtain the following exponentially 
decay estimate 



00 n 

\Ue{x)\<Y, E 

n=0 m=—n 



2(2n + 1) 



k^Rih^^\kR{] 



Mae" 



as £ — )■ +0, by using f l4.23p . fl4.24p and fl4.26p . This verifies the sharpness of 
Theorem 11.31 



Appendix 

We shall give a proof of the well-posedness of the scattering problem (11. 5p . which 
was also needed in the proof of Lemma 13.31 We could not find a convenient 
literature for the results, so for completeness we present it in this appendix. 
Our argument follows the Lax-Phillips method presented in |5]. 
Let _ 

{1,1 inM^\l], 
7,g mn\D, (4.27) 

e,r]Q + iTQ mD. 

Then the scattering problem (jl.Sp can be formulated as follows: 
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Find u G Hfg^(R ) such that u = + in M \f2 and solves the equation 



'^u:=V ■{aVu) + k^/3u = f in 



|a;|-5>oo \d\x 

where we assume supp{f) C Bfi^\Q. 

The uniqueness of the solutions to the system f l4.28p can be shown in a 
similar argument as the one used in the proof of Lemma 13.11 Next we show 
only the existence and stability estimate. 

In the following, by appropriately choosing Rq we can assume that k"^ is not 
a Dirichlet eigenvalue in Bfig^i. Let 6{x) G C°°(M^) be a cut-off function such 
that 9{x) = for |x| < Rq and 9{x) = 1 for |x| > Rq + 1. Setting 

W = u in n and W = u' + {1 - 9)u' in M^\n, (4.29) 

we can then verify directly that W G HI^^{MJ^) satisfies 

' = g in M^, 



|x|-5>oo ' ' \d\x 

with ^? = -(A + e){eu^) + / G L\BR,^,\n). 

Next, we look for a solution to fl4.30p of the following form 

W = w-(t){w-V), (4.31) 

where is C°° cut-off function such that = 1 in B^^^ and </> = in MJ^\B ^^j^i. 
We let V G H^{Br^+i) be the solution of the system 




on SBr^+i 



(4.32) 



and w G Hl^J^^) be the solution of the system 

'(A + P)w = 5(* in 
hm \x\^^^^y-[^-^kn^=^, ^^"^^^ 



,|x|-s.oo \d\x 

where g* G L"^ {B R^^J^l\VL) shall be determined later. 
Clearly, by the classical regularity estimates we see 

V^H\Br,+^\^) and weHl^{m!') 
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By direct verification we have 

g =(A + k'^)W = Aw + + A(f){w - V) 

+ 2V0 ■ V{w -V) + (f) {A{w -V) + k\w - V)) (4.34) 
=9* + Kg*, 

where K is defined to be Kg* = A(t){w - V) + 2V0 ■ V(w - 1/). 

We can show that K is compact from L'^{BFto+i\^) to itself. We shall make 
use of the Fredholm theory to show the unique solvability of fl4.34p . It suffices 
to show the uniqueness of solution to fl4.34p . We set g = 0. By fl4.30p we have 
W = 0. Hence w = 0(w - V") in and V" = in and w = in R^\Br^+i. 
It is straightforward to verify that 

\A + k')iV-w) = in Bn,+u ^^^^^ 
V -w = on OBr^+i, 

hence V — w = Therefore u; = 0, which then implies that g* = 0. Then by 
the Fredholm theory we have a unique g* G L"^ {B fjg^i\Q) to fl4.34p such that 

\\g*\\L^BR^+,\n) < C\\g\\L2(Bno + i\n) < C + \\f\\LHBR^\n)^ ■ 

s Finally, by the classical theory on elliptic equations one can show that 

M\m{Br,^^+,\n) < C {\\f\\LHBn,\n) + \W'\\m{Bno+i\n)) ■ 
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